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ABSTRACT 



We study the correlations between timing and X-ray spectral prop- 
erties in the low mass X-ray binary 4U 0614+09 using a large (265-ks) 
data set obtained with the Rossi X-ray Timing Explorer. We find strong 
quasi-periodic oscillations (QPOs) of the X-ray flux, like the kilohertz 
QPOs in many other X-ray binaries with accreting neutron stars, with 
frequencies ranging from 1329 Hz down to 418 Hz and, perhaps, as low 
^3 ■ as 153 Hz. We report the highest frequency QPO yet from any low 

mass X-ray binary at 1329 ± 4 Hz, which has implications for neutron 
star structure. This QPO has a 3.5 a single-trial significance, for an esti- 
mated 40 trials the significance is 2.4 a. Besides the kilohertz QPOs, the 
Fourier power spectra show four additional components: high frequency 
noise (HFN), described by a broken power-law with a break frequency 
between 0.7 and 45 Hz, very low frequency noise (VLFN), which is fitted 
as a power-law below 1 Hz, and two broad Lorentzians with centroid 
frequencies varying from 6 to 38 Hz and 97 to 158 Hz, respectively. We 
find strong correlations between the frequencies of the kilohertz QPOs, 
the frequency of the 6—38 Hz broad Lorentzian, the break frequency of 
the HFN, the strength of both the HFN and the VLFN and the posi- 
tion of the source in the hard X-ray color vs. intensity diagram. The 
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frequency of the 97—158 Hz Lorentzian does not correlate with these 
parameters. We also find that the relation between power density and 
break frequency of the HFN is similar to that established for black hole 
candidates in the low state. We suggest that the changing mass accretion 
rate is responsible for the correlated changes in all these parameters. 

Subject headings: accretion, accretion disks - - black holes - - stars: in- 
dividual (4U 0614+09) - - stars: neutron - - X-rays: stars 
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1. Introduction 

Accretion in neutron star low mass X-ray 
binaries can be studied through the spec- 
tral and timing properties of the associated 
X-ray emission. Some of these systems ex- 
hibit quasi-periodic oscillations with frequen- 
cies ranging from a few hundred Hz to more 
than 1000 Hz (kilohertz QPOs). These oscil- 
lations have been observed with Rossi X-ray 
Timing Explorer (RXTE; for reviews and ref- 
erences see van der Klis 1998, 1999; Swank 
1998 and [http://www.astro.uva.nl/ ~ecford/ 
qpos.html). The high frequency of these sig- 
nals nearly certainly reflects the short dynam- 
ical time scale in the region near the com- 
pact object where they originate. The timing 
properties at low frequencies [y < 100 Hz) as 
well as the spectral properties are the basis 
of a classification of the systems as Z oy Atoll 
sources (Hasinger & van der Klis 1989). The 
low frequency part of the Fourier power spec- 
tra of the atoll sources is dominated by two 
features (first identified with the EXOSAT 
satellite): a power-law red noise component 
in the lowest range of the spectrum {y < 1 
Hz) called very low frequency noise (VLFN), 
and a high frequency noise (HFN) component 
which is flat at low frequency and breaks to 
follow a power-law with an index of about 1 
at higher frequency. In some cases, a low fre- 
quency Lorentzian is also seen (e.g. Wijnands 
& van der Klis 1999a). The HFN is very simi- 
lar in shape to the low state noise of the black 
hole candidates (BHCs; van der Klis 1994a; 
Yoshida et al. 1993). Also the correlation be- 
tween the break frequency of the HFN and 
the centroid frequency of the low frequency 
Lorentzian is the same for black hole candi- 
dates and neutron stars (Wijnands & van der 
Klis 1999a). 

These timing properties are connected with 



the properties of the X-ray energy spectrum. 
This is seen very clearly in the Z sources both 
for the lower frequency features (e.g., van der 
Klis et al. 1985; Hasinger & van der Klis 1989; 
Wijnands et al. 1997a; Kuulkers et al. 1997; 
Kuulkers et al. 1996) and for the kilohertz 
QPOs (van der Klis et al. 1996; Wijnands et 
al. 1997b; Jonker et al. 1998). Similar corre- 
lations exist for atoll sources as well, both for 
the lower frequency features (e.g., Hasinger 
& van der Klis 1989; Prins & van der Klis 
1997; Mendez et al. 1997) and for the kilo- 
hertz QPOs (e.g., Ford et al. 1997b; Kaaret et 
al. 1998; Mendez et al. 1999). Also connec- 
tions have been made between the frequency 
of the kilohertz QPOs and the frequencies of 
the lower frequency features (van der Klis et 
al. 1996; Ford & van der Klis 1998; Wijnands 
& van der Klis 1998). The most obvious rea- 
son for all these correlations is that the chang- 
ing mass accretion rate causes changes in the 
timing and spectral properties (Hasinger & 
van der Klis 1989). 

4U 0614+09 has been studied previously 
with EXOSAT (Singh & Apparao 1994; Bar- 
ret & Grindlay 1995). The EXOSAT data 
showed an anti-correlation of high and low 
energy emission (Barret & Grindlay 1995), 
and observation with the Burst and Tran- 
sient Source Experiment (BATSE) and the 
RXTE all-sky monitor showed that this anti- 
correlation extends to about 100 keV (Ford 
et al. 1996). X-ray bursts from 4U 0614+09 
were identified by OSO-8 (Swank et al. 1978) 
and WATCH (Brandt et al. 1992; Brandt 
1994). The only previously published study 
investigating the broadband (~0. 03-1200 Hz) 
X-ray power spectra is the Mendez et al. 
(1997) study of 14 ks of RXTE data, con- 
cluding that this source can be classified as 
an atoll source. Strong kilohertz QPOs were 
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discovered in this source with RXTE (Ford et 
al. 1997a). The frequency of these QPOs cor- 
related with count rate over short time scales, 
but on long time scales this correlation broke 
down (Ford et al. 1997a; Mendez et al. 1997). 
A more robust correlation was shown to ex- 
ist between the frequency and either the flux 
of a spectral component, the X-ray spectral 
shape, or a hardness ratio (Ford et al. 1997b; 
Kaaret et al. 1998; Mendez et al. 1997). The 
physical implication is that both these latter 
parameters and the frequency of the kilohertz 
QPOs are determined by a changing mass ac- 
cretion rate. The 265 ks of RXTE data we use 
here includes 75 ks previously used by Ford 
et al. (1997a;1997b) and Kaaret et al. (1998) 
and 14 ks by Mendez et al. (1997). 

In this paper we show a connection be- 
tween the changes in spectral shape and both 
the kilohertz QPOs and the lower frequency 
features in a large RXTE data set for 4U 
0614+09. In §2 we describe our analysis of 
the spectral and timing data. In §3 we present 
our results and in §4 we compare our results 
to other atoll sources and discuss the connec- 
tions to physical models. 

2. Observations and Data Analysis 

In this analysis we use data from RXTE's 
proportional counter array (PCA; for more in- 
strument information see Zhang et al. 1993). 
The total observing time is about 265 ks, ap- 
proximately 75% of the total RXTE observing 
time on 4U 0614+09 to date. Data sets are 
split into near continuous time intervals of ap- 
proximately 2500 s which we shall call obser- 
vations. These observations are listed in Ta- 
ble 1. We use the 16-s time-resolution Stan- 
dard 2 mode to calculate the colors (as defined 
below) and the 122— /j,s time-resolution Event 
and 0.95— /is time-resolution Good Xenon modes 



for Fourier analysis. No X-ray bursts were 
observed in this data set. 

We calculate a hard color, defined as the 
count rate in the energy band 9.7—16.0 keV 
divided by the rate in the energy band 6.0—9.7 
keV, and a soft color, defined as the count 
rate in the energy band 3.5—6.0 keV divided 
by the rate in the energy band 2.0—3.5 keV. 
We also calculate the intensity, the count 
rate over five detectors in the energy band 
2.0—16.0 keV. To obtain the count rates in 
these exact energy ranges we interpolate lin- 
early between count rates in the PCA chan- 
nels. We calculate the colors and intensity 
first for each time interval of 16 s. We sub- 
tract the background contribution in each 
band using the standard background model 
for the PCA version 2.1b. We then calculate 
the average and standard deviation for each 
observation among the 16 s points. In the left 
panel of Figure p] we show a diagram of hard 
color vs. intensity, in the right panel a dia- 
gram of hard color vs. soft color. Each point 
in the plot corresponds to the average of the 
16 s points in one observation. The error bars 
indicate the standard deviation from the 16 s 
intervals. In the few observations where only 
three or four of the five PCUs were on, we 
normalize the count rates to five detectors. 

Most of the data fall on one curve. How- 
ever, there is a set of points in the 0.4 — 0.5 
hardness range that appears to be shifted to 
higher count rate. These are the points from 
the 1997 observations. There is no obvious 
reason for the shift. The general trend in the 
soft color vs. hard color diagram is the same 
as in the diagram of hard color vs. intensity 
including the offset of the 1997 points. We 
use the diagram of hard color vs. intensity 
in this source instead of the color-color dia- 
gram, as was done with 4U 1608—52 (Mendez 
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et al. 1999) and 4U 1728-34 (Mendez & van 
der Klis 1999), because there is less scatter in 
count rate than in soft color. 

Most of our observations occurred in gain 
epoch 3 (between 1996 April 15 and 1999 
March 22). Four of our observations (a to- 
tal of 9.4 ks of data) took place in gain epoch 
1 (before 1996 March 21), two (5 ks) in gain 
epoch 2 (1996 March 21 - 1996 April 15) and 
two (6.9 ks) took place in gain epoch 4 (1999 
March 22 to present). In order to correct the 
colors for these data for the changes in instru- 
mental response, we calculate colors of the 
Crab nebula, which can be supposed to be 
constant in its colors, using the same energy 
bands in the various gain epochs and calcu- 
late a correction factor which we apply to the 
data in epochs 1, 2, and 4 to make the Crab 
colors match those of epoch 3 (c.f. Kuulkers 
et al. 1994). The correction factor for hard 
colors for epochs 1,2 and 4 is about 0.2 %, 
1.0 % and 0.2 % respectively. The correction 
factor for soft colors for epochs 1,2 and 4 is 
about 3.5 %, 1.0 % and 20 % respectively. Fi- 
nally, the correction factor for intensities for 
epochs 1,2 and 4 is about 1.0 %, 6.0 % and 
2.0 % respectively. We have tested the va- 
lidity of these corrections for epochs 1 and 2 
by using spectral fits of 4U 0614+09 observa- 
tions and the response matrices to calculate 
cross-epoch corrections. The correction fac- 
tors obtained in this way are within 1.5 % 
larger in hard color, within 7 % larger in soft 
color and within 10 % smaller in intensity as 
the Crab-based corrections. 

We parametrize the position of each obser- 
vation in the hard color vs. intensity diagram 
by a variable, S a , which measures the dis- 
tance along the curve traced out by the data 
(see Figure |l|). This is inspired by the work 
on Z sources where the distance along the Z 



track in the color-color diagram is parame- 
terized by S z (Hasinger et al. 1990; Hertz et 
al. 1992) and by recent results on the atoll 
sources 4U 1608-52 (Mendez et al. 1999) 
and 4U 1728-34 (Mendez & van der Klis 
1999). We note that these previous studies 
use color-color diagrams as opposed to the 
color-intensity diagram we use here. 

In calculating S & , we choose a few points 
in the hardness-intensity diagram and fit a 
quintic spline through these points and define 
S a = 1 as the point (174,0.629) and 5 a = 3 
as the point (1255,0.295) as indicated in the 
figure. We rescaled the axes so that both axes 
range over the same values. We then find the 
S a value of each observation by determining 
the closest point on the curve. The error in 
S a is obtained by quadratically adding the er- 
rors in S a that resulted from mapping the one 
sigma errors in hard color and intensity on to 
S a . 

We calculated power spectra from the high 
time resolution data. For each observation we 
make two spectra. One in which we use all the 
energy channels of the PCA, for the fit at low 
frequency [y < 400 Hz), and one in which we 
apply channel cuts to select high energies for 
the fit of the kilohertz QPOs. The QPOs are 
generally stronger at higher energies. For the 
kilohertz QPO fits we use a nominal energy 
range of 4.6—97.3 keV (unbinned PCA chan- 
nels 13—249 for gain epoch 3), for the other 
epochs we use the channel range closest to 
these energies. For gain epoch 1 we use chan- 
nels 17-249 (4.2-68.0 keV), for gain epoch 2 
channels 13—249 (3.5—75.6 keV) and for gain 
epoch 4 we use channels 11—219 (4.6—97.5 
keV). In all cases, we divide the data into 
segments of 64 or 256 s and then bin the 
data in time before Fourier transforming such 
that the Nyquist frequency is always 4096 Hz. 
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We then average the individual power spec- 
tra over each observation. Four representa- 
tive power spectra are shown in Figure 0. We 
fit each spectrum with a function that was 
found to describe all the power spectra well. 
The four main components are as follows. (I) 
The Poisson level. For the Poisson level we 
fit a value constant in frequency. This is valid 
since dead-time effects are very small in this 
weak source. This value is approximately 2 
as expected from Poisson noise with the nor- 
malization of Leahy et al. (1983). We fit the 
Poisson level from 2000 to 4000 Hz, and then 
fix it and fit the low frequency power spec- 
tra up to 400 Hz. When we fit the kilohertz 
QPOs we let the Poisson level float. (II) A 
broken power-law whose power, P, is defined 
as: 




Av a ° [y < z/ break ) 



where A, tWik, ao and a± are the normal- 
ization factor, break frequency and the index 
below and above the break. Following previ- 
ous work we will refer to this component as 
high frequency noise (HFN; Hasinger & van 
der Klis 1989). (Ill) A power-law below 1 
Hz called very low frequency noise (VLFN; 
Hasinger & van der Klis 1989). (IV) The kilo- 
hertz QPOs, described by Lorentzians with 
best-fit centroid frequencies from 153 to 1329 
Hz. In one power spectrum the centroid fre- 
quencies of the QPOs fall in the same fre- 
quency range as the low frequency noise [y < 
400 Hz), in this case we include the shape of 
the low frequency noise in our fit of the kilo- 
hertz QPOs. 

In some cases we add two additional com- 
ponents to the fit function : a Lorentzian 
with a best-fit centroid frequency between 6 
and 38 Hz (FWHM 2-19 Hz), and another 
Lorentzian with a best-fit frequency between 



97 and 158 Hz (FWHM 22-144 Hz). Other 
atoll sources show similar components in their 
power spectra, as seen with GINGA in 4U 
1608-52 (Yoshida et al. 1993) and with 
RXTE in 4U 1705-44 (Ford, van der Klis & 
Kaaret 1998), 4U 1728-34 (Ford & van der 
Klis 1998), 4U 1820-30 (Zhang et al. 1998a), 
4U 1608-52 (Mendez et al. 1999; Ford et al. 
1999) and in other binaries as well (Wijnands 
& van der Klis 1999a; Psaltis et al. 1999). 
The x 2 /dof of our fits are typically 1.0 to 1.4 
with 127 to 164 degrees of freedom. 

3. Results 

The broken power-law component of the 
power spectra is similar to that in previ- 
ous measurements in other atoll sources and 
black hole candidates. The rms fractional 
amplitude (0.1—400 Hz) ranges from 3 to 30 
%. The average value of a is close to zero 
(—0.02), and the average value of a± is close 
to unity (0.92). There is no evidence for a 
second break in this component as is seen 
in some descriptions of the power spectra of 
the BHCs and the atoll source 4U 1608-52 
(Yoshida et al. 1993). Above S a values of 
2, the broken power-law has a peaked shape 
(Hasinger & van der Klis 1989): ao is between 
—0.2 and —1.6 and «i is between 1 and 2. 
These power spectra could also be fitted with 
a broken power-law with ao ~ and a\ ~ 1 
if we add an extra Lorentzian with a centroid 
frequency near z/break- The fits reported here, 
however, do not include this Lorentzian as it 
is less than 3a significant in most cases. 

Historically other functions have been used 
to describe this HFN component: a cut-off 
power-law described by P oc v~ a exp~^ c 
(Hasinger & van der Klis 1989), or Lorentzians 
with a central frequency of zero Hz (Grove 
et al. 1994; Olive et al. 1998). We have 
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tried such functions, but the F-test statis- 
tics show no preference for one of the func- 
tions. We employ the broken power-law be- 
cause it uses fewer free parameters (with the 
other functions we often need to add an extra 
Lorentzian with a frequency near z/ break ). Fits 
of the correlations (see below) also have better 
X 2 /dof with the broken power-law function. 

The VLFN is present in some of the inter- 
vals of March, November and December 1998 
(see Table 2 and Fig. |2|) with rms fractional 
amplitude (0.01—1 Hz) between 0.5 and 3.2 
%, when S a > 2.49 (see Table 2 and Fig. g). 
This state (extremely soft with VLFN) is of- 
ten referred to as the banana state (Hasinger 
& van der Klis 1989). In all but one of the 
observations where the VLFN was present, it 
is accompanied by a significant HFN compo- 
nent with rms fractional amplitude between 
3.1 and 13.6 % (see Table 2). 

The kilohertz QPO parameters (fitted in 
the energy range 4.6-97.3 keV) are listed in 
Table 1. The kilohertz QPOs have similar 
properties to those in other low mass X-ray 
binaries. The lower frequency QPO ranges 
from 824 Hz down to 418 Hz and, perhaps, 
as low as 153 Hz. The upper frequency QPO 
ranges from 1329 down to 449 Hz. Note that 
these are fits of the power spectra above 4.6 
keV. We include the kilohertz QPOs when 
the significance based on the error in the in- 
tegrated power is above 3.0 o. Of the kilo- 
hertz QPOs reported here, 69 % have a sig- 
nificance of 3.5 to 8.9 a and 31 % have a sig- 
nificance of 3.0 to 3.5 a. The QPO with the 
highest observed centroid frequency, 1329 ± 
4 Hz is 3.5 er significant. This significance is 
for a single trial (van der Klis 1989). We es- 
timate that the number of trials is about 40 
based on the width of the detection interval 
(~ 1000 Hz) and the width of the QPO (~ 



25 Hz) at 1329 Hz. Including these trials re- 
duces the significance to 2.4 a. We note that 
the QPO frequency of 1329 Hz deviates from 
the linear trend seen for the other kilohertz 
QPOs (Fig. |j). The break frequency in the 
same observation also deviates from the lin- 
ear trend seen for the break frequencies at S & 
below 2.3 (Fig. The frequency difference 
between the two kilohertz QPOs is constant 
over a wide range in frequency in 4U 0614+09 
and has a weighted average value of 311.8 ± 
1.8 Hz (Fig. f§). At least one QPO is present 
over basically the entire range of S a (1.25 to 
2.89) and in 74 of the 94 power spectra. At 
5* a > 2.49 a QPO is present in 4 of 9 observa- 
tions with rms fractional amplitudes ranging 
from 5.1 to 11.1 %. There is generally only 
one QPO above 5* a of about 2.3 and below 5* a 
of about 1.7. There is one exception at S a 
= 1.25, where there is a QPO present at 153 
± 6 Hz. However this QPO can be identified 
either as the lower kilohertz QPO or one of 
the 97—158 Hz Lorentzians (the circle point 
in Fig. |4] and fty on the basis of the corre- 
lations of centroid frequency, rms fractional 
amplitude or FWHM with S a . In most cases 
where single QPOs are observed, they can be 
identified as the upper (or, in one case lower) 
one on the basis of the correlations of cen- 
troid frequency with S a and/or hard color. 
The four highest frequency features are sin- 
gle QPOs but are likely the higher frequency 
of the two QPOs based on the S & correla- 
tion. There is, however, some ambiguity in 
this identification. 

In 31 of the total of 94 observations we 
have found it necessary to include a low fre- 
quency Lorentzian at 6 to 38 Hz. We include 
it in the fit when the significance based on 
the error in its integrated power is above 3.0 
o. This Lorentzian is seen only at low S a 
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and increases in frequency with S & (Fig. ||). 
The rms fractional amplitude (2.4—15.9 %) 
decreases as the centroid frequency increases 
(Fig. H). In some cases the Lorentzian is 
broad, in the 1996 August 8 interval for in- 
stance it has Q(=z//FWHM) of 1.3 ± 0.5. 

The final component of our fit, a Lorentzian 
near 100 Hz, is clearly present at a high S a 
value (Fig. ^j. We include it in the fit when 
the significance based on the error in its inte- 
grated power is above 3.0 a (15 observations). 
The best-fit centroid frequencies are between 
97 and 158 Hz and the rms fractional ampli- 
tudes between 3.9 and 13.8 %. Adding it to 
the fit makes a.\ and z+ rea k of the HFN signif- 
icantly larger. 

Most of the features in the Fourier spec- 
tra discussed above are correlated smoothly 
with the changing energy spectrum. This is 
summarized in Figures ^|, [5] and |41 where we 
plot several quantities vs. S a , which measures 
the changing colors. The same quantities do 
not correlate smoothly with count rate (in the 
energy band 2.0—16.0 keV). The only excep- 
tion to these correlations is the frequency of 
the 97—158 Hz Lorentzian, which does not 
correlate with S a (Fig. §). The 97-158 Hz 
Lorentzians do not fit on an extrapolation of 
the S a vs. frequency correlation of the 6—38 
Hz Lorentzians or its harmonics, suggesting 
that the 97—158 Lorentzian has a different 
physical origin. 

Since many of the parameters of the power 
spectra fits are correlated with S a , many of 
these values are correlated with each other. 
Two such correlations are shown in Figure 
[]. Both the centroid frequency of the 6—38 
Hz Lorentzian, ^lfloi-; and the break fre- 
quency of the HFN are correlated with the 
centroid frequency of the higher frequency 
kilohertz QPO, t'kHz- Fitting a relation vlflov 



= Az/£ Hz , we obtain a = 2.46 ±0.07 and A = 
3.0(±1.4) x 10~ 6 (line in Fig. |). The fit has 
a large x 2 /dof of 11.6. The quoted errors use 
A X 2 = 1.0. 

When we plot centroid frequency of the low 
frequency Lorentzian vs. the break frequency, 
the points fall in the same region as the plot 
that Wijnands & van der Klis (1999a) estab- 
lished for the BHCs, atoll sources and the mil- 
lisecond pulsar SAX J1808.4-3569. If we es- 
timate the frequency of the lower kilohertz 
QPO by subtracting the average frequency 
difference from the upper kilohertz QPO we 
can plot this vs. the centroid frequency of 
this Lorentzian. We fit this correlation and 
find that to within the errors it is consistent 
with the correlation that Psaltis, Belloni & 
van der Klis (1999) found for neutron stars 
and black holes. 

In order to compare the HFN component 
with that seen in BHCs (Mendez & van der 
Klis 1997), we also calculated the power den- 
sity at the break of the broken power-law. 
The values from 4U 0614+09 and from the 
BHCs are plotted in Figure 0. 

4. Discussion 

4.1. Timing— Spectra Correlations and 
States 

The previous section demonstrates connec- 
tions between properties of the energy spec- 
trum and features in the power spectra of the 
low mass X-ray binary 4U 0614+09. There 
is a strong correlation between the position 
in the X-ray hard color vs. intensity dia- 
gram, the frequencies of the kilohertz QPOs, 
the frequency of a 6—38 Hz Lorentzian, the 
break frequency of the HFN and the strength 
of both the HFN and the VLFN. 

Historically the atoll sources (4U 0614+09 
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among them) are known to show two source 
states, the so-called 'banana' and 'island' states 
named after their shapes in the color-color di- 
agram (Hasinger & van der Klis 1989). In the 
banana state there is a strong VLFN compo- 
nent present in the power spectra at low fre- 
quencies and the energy spectrum is soft. In 
the island state the power spectrum is dom- 
inated by the HFN and the energy spectrum 
is harder. In the present data we observe 
a smooth transition between the two states, 
indicating that these are part of a continu- 
ous sequence (Hasinger & van der Klis 1989). 
Similar transitions are seen in RXTE data 
of other atoll sources (4U 1636—53, Prins 
k van der Klis 1997; 4U 1705-44, Ford, 
van der Klis & Kaaret 1998; 4U 1820-30, 
Zhang et al. 1998a; 4U 1735-44, Ford et al. 
1998; 4U 1608-52, Mendez et al. 1999; 4U 
1728-34, Mendez & van der Klis 1999; Aql 
X-l, Reig et al. 1999; Mendez 1999). In 
4U 0614+09 the kilohertz QPOs are present 
in the whole range of the X-ray color-color 
diagram in contrast to data sets on similar 
sources that indicate that the QPOs are ab- 
sent at the highest S a and possibly the lowest 
as well (4U 1608-52, Mendez et al. 1999; 4U 
1705-44, Ford, van der Klis & Kaaret 1998; 
4U 1636-53, Zhang et al. 1996, Wijnands et 
al. 1997c; 4U 1735-44, Wijnands et al. 1998, 
Ford et al. 1998; 4U 1820-30, Smale, Zhang 
& White 1997; KS 1731-260, Wijnands & 
van der Klis 1997d). At the highest S & (> 
2.49) QPOs are present in 4U 0614+09 with 
rms fractions of 4 to 8 % in the total en- 
ergy range (5 to 11 % for an energy range of 
4.6—97.3 keV). These detections are similar to 
previously reported upper limits for the rms 
fraction over the full energy range in the ba- 
nana state: < 2-3% (4U 1735-44, Wijnands 
et al. 1998b), < 2 - 6% (4U 1705-44; Ford, 
van der Klis & Kaaret 1998) and < 1 - 5% 



(4U 1608-52, Mendez et al. 1999). In 4U 
0614+09 we are probably not seeing the up- 
per banana state, where there is no HFN and 
the VLFN becomes very strong. 

The correlation between the kilohertz QPO 
frequencies and S a in this source is similar to 
that in the Z sources with S z (van der Klis et 
al. 1996; Wijnands et al. 1997b; Jonker et 
al. 1998). The strength and frequency of the 
signals at lower frequency also show similar 
correlations. S a is a much better indicator for 
the QPO frequencies than count rate, where a 
correlation with the QPO frequencies on long 
time scales is lacking (4U 0614+09, Ford et 
al. 1997a; Aql X-l, Zhang et al. 1998b; 
4U 1820-30, Zhang et al. 1998a, Kaaret et 
al. 1999; 4U 1608-52, Mendez et al. 1999; 
and other sources, see Mendez 1999). Mod- 
els explaining the kilohertz QPOs predict, for 
different reasons, that the centroid frequen- 
cies increase with mass accretion rate (Miller, 
Lamb & Psaltis 1998; Titarchuk, Lapidus & 
Muslinov 1998; Stella & Vietri 1998). 

The present data indicate a link between 
S a and mass accretion rate. The correla- 
tion between break frequency of the HFN and 
S a further strengthens this link. The sug- 
gestion that the break frequency tracks the 
mass accretion rate was previously noted for 
atoll sources (van der Klis 1994a; Prins et 
al. 1997). The only possible exceptions to 
a correlation between break frequency and 
inferred mass accretion rate so far are SAX 
J1808.4-3658 (Wijnands & van der Klis 1998) 
and SLX 1735-269 (Wijnands & van der Klis 
1999b) where at higher count rates the break 
frequency is smaller. The exception to the 
correlation in these sources may be related 
to the fact that they were observed at espe- 
cially low inferred mass accretion rates. The 
two spectra of SLX 1735—269 have exception- 
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ally low break frequencies of 2.3 and 0.11 Hz 
(Wijnands & van der Klis 1999b). We note 
however, that accretion rates are inferred only 
from the fluxes and the relative change in flux 
is small, about 20 % in SLX 1735-269 and 
about 40 % in SAX J1808.4-3658. The anti- 
correlation in these sources may be due to this 
small dynamic range. 

The changes in energy spectra, reflected in 
a changing S a , are likely driven by changes 
in the mass accretion rate, as previously sug- 
gested for 4U 0614+09 by the correlation be- 
tween the frequency of the kilohertz QPOs 
with the blackbody flux (Ford et al. 1997b) 
or with the index of the power-law (Kaaret 
et al. 1998). Larger values of S a correspond 
to softer spectra, higher fluxes, and larger 
inferred mass accretion rates. From color- 
color diagrams, at low inferred mass accre- 
tion rate there is an anti-correlation of flux 
and hardness, which is manifest also in spec- 
tral studies of atoll sources: e.g., 4U 1636—53 
(Breedon et al. 1986), 4U 1735-44 (Smale 
et al. 1986), 4U 1705-44 (Langmeier et al. 
1987), 4U 1608-52 (Mitsuda et al. 1989), 4U 
0614+09 (Barret & Grindlay 1995; Ford et al. 
1996), and others (van der Klis & van Paradijs 
1994). One long-standing explanation of this 
anti-correlation is in terms of thermal Comp- 
tonization models (e.g. Sunyaev & Titarchuk 
1980). In the thermal model the increasing 
soft X-ray flux cools the Comptonizing elec- 
tron plasma, leading to a decrease in its tem- 
perature and thus a decrease in hard X-ray 
flux. 

4.2. Comparison with Other Sources 

The VLFN component might correspond 
to motion of the source along the track in 
the color-color diagram during an observa- 
tion, as suggested for Z sources (van der Klis 



1991). Observations of Cyg X— 2 however in- 
dicate that the motion in the Z track is at 
most partly responsible for the VLFN (Wi- 
jnands et al. 1997a). Bildsten (1995) sug- 
gested that the VLFN is caused by unstable 
burning of nuclear fuel in a regime of time de- 
pendent helium burning. The dependence of 
the strength of the VLFN on the mass accre- 
tion could then be explained by a slower mode 
of combustion that sets in at higher mass ac- 
cretion rate. This slower mode of combustion 
leads to stronger long time scale variations. 

The 6— 38 Hz Lorentzians are similar to the 
horizontal branch oscillations (HBOs) in the 
Z sources in their rms fractions and correla- 
tions with S a and z/kHz- The HBOs sometimes 
show a saturation of the HBO frequency with 
respect to S z at a frequency of about 60 Hz 
(GX 17+2, Wijnands et al. 1997b; Cyg X- 
2, Wijnands et al. 1998a). In 4U 0614+09 
the Lorentzian frequency only goes up to 38 
Hz and we do not observe a saturation of the 
Lorentzian frequency with respect to S a (Fig. 
|5]), but we do observe a saturation of the break 
frequency of the HFN with respect to (Fig. 

D- 

The 97—158 Hz Lorentzians could be re- 
lated to the 65—317 Hz oscillations seen in 
some of the black hole candidates (GRO J1655- 
40, Remillard et al. 1999a; XTE J1550-564, 
Remillard et al. 1999b; GRS 1915+105, Mor- 
gan et al. 1997). These oscillations fall 
in the same frequency range as the 97—158 
Hz Lorentzians we observe but have weaker 
rms fractional amplitudes (~1 %). In GRS 
1915+105 (65-68 Hz; Morgan et al. 1997) 
and GRO J1655-40 (281-317 Hz; Remillard 
et al. 1999a) the oscillations were about con- 
stant in frequency. In XTE J1550-564 how- 
ever, the frequency varied considerably from 
184 to 284 Hz (Homan et al. 1999; Remillard 
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et al. 1999a). 

For what concerns the other features in the 
power spectra, the high frequency noise com- 
ponent in atoll sources is very similar to the 
low-state noise in black hole candidates (van 
der Klis 1994a,b). The BHCs show a correla- 
tion of power density at the break with break 
frequency (Belloni & Hasinger 1989). The re- 
lation between power density at the break and 
break frequency is similar in 4U 0614+09 and 
the BHCs in the low state (Fig. [F]). However 
in the region where the break frequencies of 
4U 0614+09 and the BHCs become compara- 
ble (the intermediate and very high state of 
the BHCs) 4U 0614+09 does show a larger 
power density. The HFN may result from a 
superposition of shots with the break corre- 
sponding to the shots with the longest dura- 
tions, as proposed for the BHCs (e.g. Belloni 
& Hasinger 1990). The time scale of the shots 
could be set by the lifetime of clumps (van der 
Klis 1994b) and could decrease with mass ac- 
cretion rate (Ford & van der Klis 1998). 

4.3. Correlations of Timing Features 

Stella & Vietri (1998) have proposed that 
the 6—38 Hz Lorentzian is caused by relativis- 
ts frame dragging in the inner accretion disk 
(Lense-Thirring precession). The frequency 
of the Lorentzian depends on the frequency 
of the upper kilohertz QPO, assumed to be 
fixed by the Keplerian frequency at the in- 
ner edge of the disk. Stella & Vietri (1998) 
predict a scaling index of 2, but see correc- 
tions by Morsink & Stella (1998). Our value 
of 2.46 ± 0.07 is significantly different from 
2. Also, as in 4U 1728-34 (Ford & van der 
Klis 1998), the required ratio of the moment 
of interia over the neutron star mass is un- 
physically large if the nodal precession fre- 
quency is the QPO frequency. However this 



is not a problem if the QPO frequency is twice 
the nodal precession frequency (Stella & Vi- 
etri 1998). 

Titarchuk & Osherovich (1999) associate 
the 6—38 Hz Lorentzian with radial oscilla- 
tions in a boundary layer where the disk ad- 
justs to the rotating neutron star and the 
break frequency of the HFN with the diffu- 
sion time scale in this region. The correlations 
in 4U 0614+09 of the lower kilohertz QPO 
vs. 6—38 Hz Lorentzian and HFN break fre- 
quency do not fall on the predicted curves. In 
fact the correlations with respect to the higher 
frequency kilohertz QPO are closer to those 
predicted. In addition there is a high fre- 
quency turnover of the break frequency which 
is not predicted (Fig. 

4.4. The Highest Frequency QPOs 

We report here the highest frequency QPO 
of any low mass X-ray binary to date at 1329 
± 4 Hz. Such a high frequency oscillation 
puts new constraints on the structure of neu- 
tron stars (Miller, Lamb & Cook 1998). The 
mass and radius are bounded by: 

M < 2.2(z/^/1000 Hz) -1 (l + 0.75j) M Q 

and 

R < 19.5(^/1000 Hz) -1 ^ + 0.2 j) km 

where j is the dimensionless angular momen- 
tum of the star, j = 1ixcIv svin IGM 2 (e.g. 
Miller, Lamb & Psaltis 1998). This expan- 
sion in j is valid for spin frequencies less than 
about 400 Hz (Miller, Lamb & Cook 1998). 
For I = 3xl0 45 g cm 2 (Cook, Shapiro & 
Teukolsky 1994; Miller, Lamb & Psaltis 1998) 
and v S pin — 311.8 Hz, the inferred spin fre- 
quency of the star, the kilohertz QPO at 1329 
± 4 Hz constrains M and R to M < 1.9M m 
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and R < 15.2 km. These constraints rule out 
the extremely stiff EOS. A mean-field theory 
('L' in Cook, Shapiro & Teukolsky 1994) is 
ruled out and the tensor interaction EOS ('M' 
in Cook, Shapiro & Teukolsky 1994) is ruled 
out except for a very small range of R and 
M. However we note, that this QPO has a 
single-trial significance of only 3.5 a. To put 
constraints on the softer EOS the kilohertz 
QPO must be observed at higher frequencies. 
For a slow rotating star, 2010 Hz would rule 
out the 'U'-type EOS (Wiringa, Fiks & Fab- 
rocini 1988) and 1920 Hz would rule out the 
commonly used 'FPS' EOS (Cook, Shapiro & 
Teukolsky 1994). 

We note that there is no obvious saturation 
of the kilohertz QPO frequency with respect 
to S a , our chosen indicator of mass accre- 
tion rate (see Fig. ^). However the kilohertz 
QPO at 1329 Hz does deviate from the lin- 
ear trend seen for the other kilohertz QPOs. 
Such a saturation where the frequency be- 
comes independent of the chosen mass accre- 
tion rate indicator would be an indication for 
the marginally stable orbit predicted by gen- 
eral relativity (Kluzniak, Michelson & Wag- 
oner 1990; Kaaret, Ford & Chen 1997; Miller, 
Lamb & Psaltis 1998). Such a saturation has 
been reported in 4U 1820—30 (Zhang et al. 
1998a; Kaaret et al. 1999). 

4.5. QPO Frequency Difference 

We found that, Az/, the frequency differ- 
ence between the upper and the lower kilo- 
hertz QPO peak is constant in 4U 0614+09, 
confirming the results of Ford et al. (1997a). 
This matches the prediction of the beat fre- 
quency mechanism (e.g. Miller, Lamb & 
Psaltis 1998) in which Av is equal to the spin 
frequency and therefore constant. However 
observations of five sources have shown a sig- 



nificant decrease in Av as the frequency of the 
kilohertz QPOs increases (Sco X-l, van der 
Klis 1997; 4U 1608-52, Mendez et al. 1998a,b; 
4U 1735-44, Ford et al. 1998; 4U 1728-34, 
Mendez & van der Klis 1999; 4U 1702-42, 
Markwardt et al. 1999a,b). To examine the 
trend of Av with increasing QPO frequency 
in 4U 0614+09 we plot Av vs. the frequency 
of the lower frequency kilohertz QPO for all 
these sources in Figure |8|. We find that at 
low QPO frequency (< 750 Hz) Av in 4U 
0614+09 is consistent with Sco X-l and 4U 
1608-52 and at high QPO frequency (750-820 
Hz) inconsistent (2.8-8.0 o different). How- 
ever, we note that one of the points from 4U 
0614+09 itself deviates 11.0 a from the trend 
seen in 4U 0614+09. 
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Table 1: 4U 0614+09 RXTE Observations & kilohertz QPO fits 



T re3 Rate Higher kHz QPO Lower kHz QPO 

frequency FWHM rma frequency FWHM 

(CO (=W0 (Hz) (Hz) (%) (Hz) (Hz) 



1.14±0.1C 

i.ie±o.o9 

L.23±0.0S 
l.24±0.08 
l.24±0.09 
l.25±0.08 
L.26±0.09 
L.28±0.08 
l.30±0.09 
l.32±0.07 
l.32±0.07 
l.36±0.07 
l.3S±0.07 
l.39±0.07 
l.40±0.07 
l.44±0.07 
l.44±0.07 
l.44±0.06 
l.44±0.07 
l.44±0.08 
l.49±0.08 
l.50±0.06 
l.51±0.06 
l.51±0.06 
l.51±0.07 
l.53±0.0S 
l.55±0.0S 
L.56±Q.QS 

l.57±0.07 
l.58±0.08 
l.58±0.07 
l.58±0.07 
l.58±0.06 
l.60±0.0S 

i.6i±o.os 

L.61±0.07 



2fi Feb 
2 6 Feb 
16 Mar 



23 De< 
22 Apr 
22 Ap: 
18 0f 



10 Jan 
10 Jan 
13 Mar 

6 Sep 



1Q96 08:64 
1996 10:34 
1998 11:45 

1996 03:32 
1998 10:02 

1996 01:60 
1998 11:43 

1996 00:11 
1998 10:23 

1997 13:28 
1996 22:31 

1996 20:55 
1996 19:19 

1998 21:33 
1996 19:21 2 

1996 22:41 

1997 18:21 
1999 23:47 

1996 20:58 2 

1998 11:43 
1999 22:06 

1997 10:49 2 

1997 10:16 

1998 10:25 
997 22:55 



> 19' 



30 Dec 
13 Apr 
15 Mai 
18 Mai 

2 6 May 
6 Sep 
3 Sep 



1998 22:07^ 

1996 11:3a 1 
1998 11:44 
1998 13:35 

1999 OS:0S 2 

1997 17:54 
1997 21:12 
597 19:28 



i3 



S55.6±14.9 
S85.5±32.1 
547.3±11.9 

569.2±25.S 
S32.1±14.0 
626.1±12.2 
S75.9±20.5 
695.1±10.1 



600.2±19.9 
657.6±3.5 

S91.4±16.3 
654.6±7.6 

sea.a±29.3 

644.5±33.7 



706.9±19.0 
605.9±17.7 
623.5±14.2 
626.8±13.3 
5S1.5±16.2 



103.6±50.6 
293±139 
72.0±37.8 

196.6±76.9 
126.2±77.9 

S4.6±31.2 
131.2±44.0 

84.8±26.6 



211.0±93.4 
19.9±8.1 

133.1±57.9 
37.7±28.5 
252±125 

146.7±49.0 
235±107 



135.5±60.9 
136.2±54.S 
167.7±55.S 
178.6±46.3 
115.1±45.2 



16.7±3.1 
24.9±5.3 
13.1±2.6 

19.3±3.S 
15.4±3.5 
18.0±2.9 
1S.5±2.S 
19.0±3.0 



21.4±4.2 
12.7±1.8 
16.5±2.8 
12.0±2.5 
15.6±3.6 
16.0±2.4 
22.5±4.9 



19.1±3.8 
17.6±2.8 
17.8±2.5 
19.3±2.3 
1S.1±3.1 



l.61±0.06 
L.63±0.07 
l.65±0.0S 
l.67±0.04 
l.67±0.04 
l.68±0.0£ 
l.71±0.0S 
l.71±0.04 
l.72±0.06 
l.77±0.04 
l.79±0.04 
l.79±0.05 
l.81±0.0S 
l.83±0.04 
l.83±0.04 
l.83±0.0S 
l.84±0.06 
l.84±0.0S 
l.85±0.0S 
l.86±0.04 
l.86±0.04 
l.87±0.04 
l.88±0.04 
l.90±0.05 
l.94±0.04 
l.95±0.04 
l.98±0.04 
l.99±0.04 
2.00±0.04 
2.01±0.0S 
2.01±0.0S 

a.oi±o.o4 

2.04±0.05 
2.04±0.04 
2.05±0.04 
2.06±0.03 
2.09±0.04 
2.13±0.03 
2.14±0.04 
2.14±0.03 
2.16±0.04 
2.1S±0.04 
2.20±0.03 
2.21±0.02 
2.22±0.02 
2.25±0.02 
2.25±0.02 
2.25±0.02 
2.49±0.02 
2.58±0.02 
2.59±0.03 



2.78±0.06 
2.79±0.06 
2.79±0.04 
2.89±0.06 
2.91±0.06 



12 Nov 
23 May 
16 Mar 
26 Jan 

15 Aug 

16 May 
14 Aug 
14 Aug 
19 Mar 
25 Apr 
16 Mar 
24 Dec 
24 Ap: 

6 Sep 
5 Sep 
31 De- 



1998 13:46 

1999 00:2s 1 
1997 22:25 
1997 13:39 
1997 00:06 
1997 05:23 
1997 20:51 

1997 22:32 

1998 10:04 
1996 00:11 



1996 09:11 
1998 21:35 
1996 22:36 
1997 00:16 
1997 21:05 
1998 22:07 
1996 03:31 
1996 01:47 



5 Srp 1997 22:39 
25 Apr 1996 03:23 
7 May 1997 10:14 
24 Apr 1996 20:59 
24 Apr 1996 19:48 

7 Aug 1996 01:60 
24 Apr 1996 14:58 
24 Apr 1996 18:19 

7 Aug 1996 00:10 

24 Apr 1996 19:23 
24 Apr 1996 13:21 2 

25 Apr 1996 04:59 
14 Mar 1998 11:43 2 

24 Apr 1996 16:36 
10 Mar 1998 13:24 

6 Aug 1996 22:30 
29 Dec 1998 22:07 



9 Mar 1998 10:02 
7 Mar 1998 10:01 
1 Nov 1998 15:24 
10 Mar 1998 10:02 



9 Nov 1908 1.2:2-! 
30 Oct 1998 13:67 

8 Mar 1998 10:20 
6 Nov 1998 15:23 2 

25 Dec 1998 21:37 

2 Nov 1998 17:00 4 
6 Nov 1998 1S:23 4 
8 Nov 1998 13:66 4 
31 Oct 1998 19:03 5 

3 Nov 1998 10:61 5 

4 Nov 1998 18:43 5 
27 Dec 1998 21:35 2 

6 Mar 1998 12:00 
28 Dec 1998 22:06 

26 Dec 1998 21:36 

7 Nov 1998 15:23 



702.7±11.9 
668.4±15.7 
673.7±21.1 
629.2±13.7 
644.4±8.6 
G52.5±8.3 
663.9±8.5 
67S.0±6.8 
799.1±15.1 
767.2±9.6 
720.8±10.7 
754.0±11.9 
809.4±4.2 
751.8±7.7 
750.4±10.6 
831.0±12.3 
S51.3±23.4 
803.6±5.1 

S4S.7±6.4 
777.3±11.4 
818.9±9.3 
736.3±14.1 
845.7±10.0 
B53.1±10.7 
909.7±13.1 
SSl.SiS.l 
894.3±3.6 
973.8±13.8 
871.3±7.2 
S94.7±6.5 
911.9±5.2 
1024. 7±6. 3 
934.8±7.4 
1024.7±16.2 
987.0±16.2 
973.7±20.6 
1034.4±6.1 
1106.0±22.7 

1066.6±10.4 
98S.S±5.8 
1161. 8±4. 6 
1136. 6±6. 2 
1175. 3±6. 5 
1038.6±7.1 

1016.8±11.1 
1146. 0±3. 4 
1221. 4±6. 2 
1217.2±3.1 

1272.6±14.6 



72.6±46.7 
130.2±73.8 

224±103 
156.0±73.1 

94.6±27.0 
123.5±37.C 
133.3±28.9 

92.0±29.4 
102.1±36.1 
134.7±33.4 
168.8±38.8 
152.8±73.0 

39.3±14.5 
106.1±22.9 
155.6±36.0 
167.2±61.4 
1S5.7±77.S 

66.0±13.5 

81.8±22.6 
129.4±27.3 
145.2±32.7 
116.6±49.6 
115.1±29.0 
101.4±37.8 
96.1±82.6 
66.4±23.7 
26.6±12.5 
132.4±37.4 
39.6±11.0 
46.6±19.9 
66.8±17.4 
41.1±30.8 
67.1±18.7 
109.2±54.0 
81.0±31.6 
341±109 
34.1±32.1 
230.6±96.6 

71.4±29.1 
39.4±21.3 
66.7±11.S 
46.2±15.5 
S3.7±27.2 

21.5±8.1 
64.6±39.4 
28.4±15.2 
62.2±13.3 

34.1±9.5 
89.6±30.5 



12.7±2. 
18.3±4 
18.3±3 
16.3±2 
13.6±1 
16.3±2 
17.2±1 
14.4±1 
16.2±2 
20.1±2 
18.9±2 
1S.1±3 
12.1±1 
16.6±1 
17.6±1 
20.6±3 
18.6±3 
15.3±1 

15.3±1 
15. Oil 
20.6±2 
13.9±2 
17. Oil 
14.8±2 
14.9±5 
13.6±1 
10.6±1 
16.4±2 
12.4±1 
12.9±1 
12.4±1 
11.0±2 
12.6±1 
13.7±2 
11.1±1 
21.2±4 
9.6±2. 
17.0±3 

10.1±1 
7.9±1. 
13.9±1 
10.9±1.E 
13.1±1 
6.7±1,5 
9.0±1.1 
9.1±1.7 
ll.lil.o 
8.4±0.7 
6.6±0.9 



523.4±5.9 
473.1±10.9 
565.3±4.4 



5S1.2±3.8 
585.8±3.6 
561.3±2.5 
570.7±3.2 
623.8±4.2 
5S2.5±3.5 
581.7±1.5 
5S0.4±3.6 
695.8±4.5 
609.4±3.1 
707.9±3.7 
679.6±1.6 
636.6±2.4 
726.7±2.9 
762.5±7.6 



27.3±22.6 
43.9±34.6 
20.4±0.3 



10.7±14.7 
18.7±6.3 
13.0±9.0 
12.3±6.9 

42.4±11.3 

1S.1±11.6 
9.0±4.2 

35.2±12.6 

44.6±12.3 
37.1±8.4 
42.6±9.7 
26.2±5.7 
31.4±6.3 
39.5±9.3 

S3.5±23.S 
2S.S±4.9 

64.2±10.6 
45.4±4.S 

3S.1±37.4 
204.2±66.4 



7.3±1.4 
7.1±1.3 
12.6±1,3 
9.8±1.4 
8.6±1.1 
10.4±1 
14.2±1 
13.9±1 
14.8±1 
12.7±1 
11.1±Q,9 
14.8±1.4 
14.0±1 
13.0±1 
14.2±1 
13.4±0 
6.3±1, 



Note. - Liated from top to bottom in order of increasing S a (position in the hard color va. intensity diagram) are, the start time (in Universal Time, 
Coordinated), T obg (the total useful observing time), T re3 (the raw time resolution of the PCA mode we use for making the Fourier power spectra), the 
background subtracted count rate over five detectors in the energy band 2.0-16.0 keV and the parameters of the kHz QPOs fitted from energies 4.6 to 97.3 
keV. 

1 in this interval only 3 PCUs were on. 

2 in these intervals only 4 PCUs were on. 

3 > 4 > 5 these intervals were added up to improve the statistics. 



Fig. -n.- 
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Table 2: 4U 0614+09 low frequency power spectra fits 









HFN 




VLFN 


Low Fr 


quency Lor 


entzian 


High Fr 


quency Lor 


ntzian 




^break 


rms 


a 


a 1 


rms a 


frequency 


FWHM 


rms 


frequency 


FWHM 


rma 




( Hl ) 


(%) 






(%) 


(Hz) 




(%) 


(Hz) 


(Hz) 


(%) 


1.14±0.1O 


0.7±0.1 


29.S±2.3 


Q.QS±Q.QS 


1.0Q±Q.Q2 




6.3±0.5 


5.9±1.2 


11.2±1.S 








1.16±C.C9 


0.7±0.1 


26.3±2.6 


-0.06±Q.O8 


1.02±0.03 




6.7±0.6 


7.6±1.6 


15.9±2.2 








1.23±0.0S 


2.1±0.3 


26.4±2.4 


0.02±0.05 


0.99±0.04 




ll.Sil.l 


11.6±3.C 


12.1±2.0 








1.24±0.0S 


l.S±0.2 


29.2±2.0 


Q.1Q±Q.QS 


0.97±0.03 




12.2±1.4 


9.3±2.9 


9.9±1.S 








1.24±C.C9 


2.2±0.2 


28.7±2.1 


Q.Q1±Q.Q4 


0.97±C.C3 




12.7±1.4 


12.5±3.S 


10.9±2.0 








1.25±0.0S 


l.S±0.2 


2S.3±1.7 


-0.03±0.05 


0.96±0.03 




11.2±1.1 


12.1±2.Q 


13.2±1.7 








1.26±C.C9 


2.5±0.3 


30.2±2.8 


-0.01±0.04 


0.81±0.02 




13.1±C.S 


7.0±4.3 


S.9±2.4 








1.28±C.C8 


1.7±0.1 


29.2±1.3 


-0.08±0.05 


0.9Q±0.03 




13.S±1.3 


11.9±2.Q 


10.9±1.S 








1.3C±C.C9 


3.2±0.4 


2S.9±1.3 


O.OS±0.04 


O.S3±Q.Q3 
















1.32±C.C7 


2.0±0.1 


2S.0±1.2 


0.02±0.03 


0.88±C.C2 




12.4±0.S 


11.9±2.4 


11.1±1.2 








1.32±C.C7 


3.4±0.5 


31.0±O.Q 


-0.01±0.05 


O.SQ±Q.Q3 
















1.36±C.C7 


3.0±0.3 


25.4±2.6 


-0.10±0.07 


o.93±o.oe 




1S.3±2.S 


17.3±11.4 


12.4±e.l 








1.38±C.C7 


5.4±0.4 


30.4±1.2 


Q.Q4±Q.Q3 


0.86±0.03 
















1.39±C.C7 


6.3±1.1 


29.4±1.4 


0.02±0.04 


0.93±0.04 
















1.40±C.C7 


2.1±0.3 


26.2±1.1 


-0.09±Q.OS 


0.83±C.C3 




14.S±1.3 


1S.9±3.6 


13.2±1.4 








1.44±0.07 


4.7±0.3 


30.2±0.7 


-0.02±0.04 


0.86±0.03 
















1.44±C.C7 


e.i±o,7 


27.6±C9 


0.06±0.03 


O.SQ±Q.Q3 




17.1±1.S 


10,S±4.S 


7.2±2.2 








1.44±C.C6 


4.1±0.2 


25.5±1.C 


Q.QQ±Q.Q2 


0.96±0.02 




16.7±0.6 


9.2±3.1 


S.6±1.2 








1.44±0.07 


e.5±o.e 


29.2±1.1 




O.S7±Q.Q3 
















1.44±0.0S 


4.6±0.3 


29.1±C.7 


Q.Q1±Q.Q4 


O.SB±0.02 




1S.S±0.4 


i.eii.i 


4.0±0.8 






- 


1.49±C.CS 


17.2±1.4 


27.5±1.6 


Q.Q7±Q.Q4 


1.00±0.05 
















1.50±0.06 


10.6±1.7 


27.2±1.1 


-o.o7±o.oe 


o.97±o.oe 
















l.siic.ce 


7.7±1.3 


29. Oil. 1 


0.04±0.05 


O.QliO.OS 
















i.5i±c.ce 


6.1±0.5 


25.6±0.Q 


0.06±0.03 


0.92±0.03 




22.6±1.2 


12.4±3.6 


S.6±1.2 








1.51±C.C7 


1S.S±2.3 


27.4±Q.S 


0.1S±0.04 


0.92±C.C6 
















1.53±C.C5 


4.9±0.4 


2S.S±Q-7 


Q.04±Q.O3 


0.91±0.03 




1S.3±0.9 


10.0±4.1 


s.oii.e 








1.55±0.0S 


e.7±i.i 


2S.S±l-2 


-o.oe±o.oe 


0.94±0.05 




21.3±C.3 


1.9±1.S 


S.l±l.l 








1.56±C.C5 


4.6±0.3 


24.2±0.7 


O.OS±0.02 


0.90±O.Q2 




1S.6±0.7 


12.8±2.6 


9.3±0.8 








1.56±C.C5 


5.7±0.4 


25.7±0.6 


0. 01±0. 02 


0.89±0.02 




20.7±0.6 


7.4±1.7 


6.7±0.e 








1.57±C.C7 


12.4±1.4 


25.3±0.9 


0.04±0.04 


0.99±C.C5 
















1.5S±C.CS 


7.8±0.9 


2S.6±1.3 


Q.Q2±0.0S 


0.96±0.05 
















1.58±C.C7 


16.4±2.C 


27.5±1.7 


Q.QS±0.04 


0.86±0.05 
















1.58±C.C7 


1S.S±1.8 


24.9±1.3 


0.06±0.04 


1.01±Q.Q6 
















l.ssic.ce 


5.6±0.5 


27.1±0.S 


0.03±0.04 


0.83±0.03 




22.9±1.1 


11.9±3.Q 


7.S±1.2 








i.eoio.os 


5.8±0.4 


24.e±0.6 


0.06±0.03 


0.88±C.C2 




20.3±1.0 


15.0±3.4 


Q.3±l.l 








1.61±0.05 


6.1±0.3 


25.e±0.6 


-O.C2±0.03 


0.8S±0.02 




24.2±C.7 


8.8±2.0 


7.3±0.7 








i.ei±o.o7 


S.0±1.1 


27.6±C.7 


-Q.02±0.04 


0.9Q±Q,04 
















1.61±0.06 


12.1±1.0 


25.4±0.6 


0.06±0.03 


0.89±0.04 
















1.63±0.07 


1C.8±1.C 


26.3±1.Q 


0.06±0.04 


0.94±Q.Q4 




26.4±0.4 


2.6±1.6 


4.7±0.9 








1.65±0.05 


6.8±0.5 


24.1±Q.Q 


Q.QQ±0.03 


0.87±C.C2 




26.6±0.9 


12.1±3.4 


7.7±0.9 








1.67±0.04 


7.5±0.6 


25.4±0.S 


0.06±0.03 


0.87±0.02 




26.6±0.7 


5.8±3.0 


S.l±1.0 








1.67±0.04 


9.2±0.4 


24.6±C.6 


0.06±0.02 


0.93±0.02 




27.0±0.7 


4.3±1.9 


4.6±0.7 








i.esio.os 


8.4±0.7 


24.7±0.6 


Q.Q4±0.03 


0.8S±0.03 




2S.3±C.9 


8.4±3.3 


S.3±l.l 








1.71±C.C5 


1C.3±C.6 


24.2±0.4 


Q.Q3±0.02 


0.93±0.02 




27.2±0.6 


6.3±1.6 


S.6±0.6 








1.71±C.C4 


9.5±0.6 


24.Q±0.4 


O.OQ±0.03 


0.94±C.C2 




29.6±0.7 


8.4±3.0 


6.9±0.8 








1.72±C.C6 


17.8±1.4 


23.0±1.4 


-C.C1±0.04 


1.00±0.05 
















1.77±C.C4 


16.4±1.3 


22.4±C8 


0.06±0.03 


1.03±0.04 










99.2±3.3 


27.0±14.3 


6.7±1.4 


1.79±C.C4 


12.8±0.8 


17.1±0.7 


0.01±0.03 


0.93±0.03 




31.6±0.7 


3.9±1.9 


3.7±0.7 








1.79±C.C5 


17.4±1.C 


21.0±0.7 


0.04±0.03 


1.02±0.04 
















l.Sl±0.05 


18.6±1.1 


22.3±0.S 


-0.03±0.03 


1.07±C.C6 










125.4±4.5 


44.0±21.1 


8.5±2.1 


l.S3±0.04 


16.6±0.9 


22.8±0.4 


-C.C5±0.04 


0.99±0.04 
















l.S3±0.04 


16.4±C.7 


21.9±C4 


-0.03±0.03 


1.01±Q.Q3 




37.6±0.6 


2.0±1.8 


2.4±0.e 






- 


l.S3±0.05 


16.3±1.6 


21.7±0.7 


-0.05±0.05 


0.87±0.05 
















l.S4±0.06 


17.4±1.S 


21.S±1.1 


-0.03±0.05 


i.oo±o.oe 
















l.S4±0.05 


18.6±1.1 


22.2±0.S 


-0.02±0.03 


0.97±CC4 
















l.S5±0.05 


1S.C±1.3 


21.8±0.S 


-C.C6±0.06 


1.03±0.07 
















l.S6±0.04 


20.6±1.3 


20.2±C7 


0.02±0.04 


1.02±0.05 
















l.S6±0.04 


16.C±C.7 


22.1±0.4 


-Q.01±0.02 


0.9S±0.03 
















l.S7±0.04 


1S.4±1.2 


21.e±0.6 


-0.02±0.04 


0.97±Q.Q4 
















l.S8±0.04 


14.6±0.8 


22.2±Q.6 


-0.04±0.03 


0.9B±0.03 
















1.90±O.C5 


2C.2±1.1 


is.e±i.o 


-0.07±0.04 


1.26±0.13 
















1.94±C.C4 


19.9±1.7 


19.8±C6 


Q.Q4±0.0S 


1.01±0.07 
















1,95±0,04 


25.0±1.5 


15.1±0.9 


-0.0e±0.05 


l.S2±Q.3S 










115.3±2.7 


25.6±12.2 


8.0±1.5 


1.98±C.C4 


2C.7±1.4 


18.4±0.6 


-o.io±o.oe 


1.07±0.08 
















1.99±0,04 




18.2±0.6 


-0.09±0,05 


















2.0C±C.C4 


26.6±1.6 


ie.o±o.s 


-C.C4±0.06 


1.3C±0.16 










15S.0±2.7 


22.3±11.1 


e.7±i.i 


2.01 i . 05 
























2.01±C,C5 






-0.12±0,08 


















2.Q 1 i . 04 


23.0^1.3 


18 6±0 S 


.0 2i0 .03 


1 .06 i . 06 
















2.04±0.05 


2S.3±3.S 


15.2±0.9 


-0.19±0.12 


1.12±0.15 
















2.04±C.C4 


23.2±1.1 


17.2±0.6 


-0.19±0.07 


1.24±0.10 
















2.0S±C.C4 


26.7±3.8 


15.4±0.S 


-0.26±0.17 


1.02±0.14 
















2.06±C.C3 


3C.S±1.2 


io.e±o.s 


-0.3e±0.12 


4.24±1.28 










119.2±ie.l 


143.7±43.7 


13.8±1.8 


2.09±C.C4 


27.9±1.C 


12.Q±Q.S 


-0.25±0.07 


2.8S±0.62 










138.0±8.3 


Q7.4±28.9 


10.9±1.3 


2.13±C.C3 


30.3±1.5 


8.2±0.5 


1.64±0.48 


5.80±1.67 










143.0±10.e 


133.5±41.3 


13.8±1.4 


2.14±C.C4 


29. £±3. 8 


13.6±0.6 


-0.2S±0.13 


1.24±0.19 
















2.14±C.C3 


31.C±2.1 


12.e±0.6 


-0.43±0.14 


1.5Q±0.20 
















2.16±C.C4 
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Fig. 1. — Hard color vs. intensity (left panel) and soft color (right panel) for 4U 0614+09. Hard 
color is defined as the ratio of count rates in the energy bands 9.7—16.0/6.0—9.7 keV. Soft color is 
defined as the ratio of count rates in the energy bands 3.5—6.0/2.0—3.5 keV. Intensity is defined 
as the count rate in the energy band 2.0—16.0 keV. Each point is an average over one observation 
(typically 2500 s). Circles in the left panel indicate values of S a along the spline curve (see text). 
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Fig. 2. — Representative power density spectra of 4U 0614+09 for four observations plotted from 
top to bottom in order of increasing S a value, S a is respectively 1.25, 1.71, 2.01 and 2.59. The 
power spectra are for the full energy band of the PCA. We subtract a constant value, somewhat 
less than the Poisson level and renormalize each spectrum to get fractional squared rms per Hz. 
The Poisson level is indicated for each spectrum with a dashed line. For clarity each spectrum 
is multiplied by a constant as indicated. Fit functions are shown which consist of a power-law 
at low frequency, a broken power-law, a low frequency Lorentzian, a Lorentzian near 100 Hz and 
kilohertz QPOs. The 25 April 1996 spectral fit does not include a ~100 Hz Lorentzian, since it 
was less than 3 a significant. 
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Fig. 3. — Parameters of the fit to the power spectra versus 5 a , the position in the hard color 
vs. intensity diagram (Fig. |IJ). The plots are rms fractional amplitude (0.1—400 Hz) of the high 
frequency noise (broken power-law), the frequency of the break in this component, and the rms 
fraction (0.01—1 Hz) of the very low frequency noise. 
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Fig. 4. — Frequency of the kilohertz QPOs vs. S a (top) and difference in frequency between double 
QPOs vs. S a (bottom). Plusses represent the higher frequency peak, squares the lower one. The 
circle marks the case where the QPO can be identified as either the lower kilohertz QPO or the 
97—158 Hz Lorentzian (see text and Fig. ||). 
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Fig. 5. — Centroid frequency and rms fractional amplitude of the 6—38 Hz (top) and 97—158 Hz 
(bottom) Lorentzian vs. S & . The circle marks the case where the QPO can be identified as either 
the lower kilohertz QPO or the 97—158 Hz Lorentzian (see text and Fig. [|). 
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Fig. 6. — Frequency of the Lorentzian at low frequencies (top) and the break frequency of the bro- 
ken power-law function (bottom) versus the frequency of the upper kilohertz QPO. The function 
shown is a power-law with index 2.46. 
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Fig. 7. — Power density of the broken power-law component at the break vs. frequency of the 
break. The filled points and the marked region are from five different black hole candidates 
(Mendez & van der Klis 1997 and references therein). The filled circles and the marked region 
are observations in the low state and the filled triangles are from the intermediate and very high 
states. The open squares are the points from 4U 0614+09. 
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Fig. 8. — The difference in frequency between double kilohertz QPOs vs. the frequency of the 
lower kilohertz QPOs for four different sources. The open triangles are the points from Sco X-l 
(van der Klis et al. 1997), the open circles are the points from 4U 1608-52 (Mendez et al. 1998), 
the open stars are the points from 4U 1735-44 (Ford et al. 1998; Wijnands et al. 1998b), the 
open diamonds are the points from 4U 1728-34 (Mendez & van der Klis 1999), the open squares 
are the points from 4U 1702-42 (Markwardt et al. 1999a,b) and the filled squares are the points 
from 4U 0614+09 (this paper). 
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